Microscopic description of d-wave superconductivity by Van Hove nesting in the 

Hubbard model 
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We devise a computational approach to the Hubbard model that captures the strong coupling 
dynamics arising when the Fermi level is at a Van Hove singularity in the density of states. We rely 
on an approximate degeneracy among the many-body states accounting for the main instabilities 
of the system (antiferromagnetism, d-wave superconductivity). The Fermi line turns out to be 
deformed in a manner consistent with the pinning of the Fermi level to the Van Hove singularity. 
For a doping rate S ~ 0.2, the ground state is characterized by d-wave symmetry, quasiparticles 
gapped only at the saddle-points of the band, and a large peak at zero momentum in the d-wave 
pairing correlations. 
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Nowadays it is widely accepted that the single-band 
Hubbard model is an appropriate starting point to de- 
scribe the electron correlations in the copper-oxide ma- 
terials jjj . In spite of the efforts to show that the model 
has a phase of <i-wave superconductivity, however, there 
is no conclusive proof at present of its existence. There 
have been different theoretical approaches to understand 
the symmetry of the order parameter, among which we 
may quote the proposals based on the effect of antiferro- 
magnetic fluctuations and the Kohn-Luttinger mech- 
anism due to Van Hove nesting of the Fermi surface . 
Numerical methods have also provided some evidence of 
enhanced pairing correlations in the 2D Hubbard model 
UK- Recent progress has been achieved by the use of 
dynamical mean-field approximations |7j] and improved 
quantum Monte Carlo methods |^||. 

We focus our discussion on the proposal of supercon- 
ducting pairing by the influence of a Van Hove singularity 
near the Fermi level pp| . Recently, some understanding 
of the problem has been attained by the use of refined 
renormalization group methods for interacting electrons 
in two dimensions |pT| , p^| . This approach suffers, though, 
from the same shortcoming of any attempt to deal with 
the origin of high-T c superconductivity. The system is 
likely to develop strong antifcrromagnetic or supercon- 
ducting correlations but, as long as the effective inter- 
action grows large as more high-energy electron modes 
are integrated out, it is not possible to discern rigorously 
the ground state of the model. A related problem con- 
cerns the fact that the superconducting correlations are 
enhanced like log 2 e, when the electron modes are in- 
tegrated out down to energy e near the Fermi surface. 
This reflects in another fashion that the model does not 
lead to a conventional scale-invariant theory at arbitrar- 
ily low energies — unless the Fermi energy is promoted to 
a dynamical quantity susceptible itself of renormalization 

In this paper we address the characterization of the 
ground state of the Hubbard model, when the Fermi level 
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FIG. 1. Points of the 8x8 lattice in momentum space. 



is at a Van Hove singularity in the density of states. For 
this purpose we will formally make use of the scaling 
properties of the model |f|, which ensure that upon in- 
tegration of the high-energy modes the electron quasi- 
particles remain in correspondence with the original one- 
particle states, up to the energy in which the effective 
interaction starts to grow large. To elucidate the be- 
havior of the system at lower energies we undertake a 
numerical diagonalization of the Hubbard hamiltonian, 
in the Hilbert space truncated to include the states left 
after integration of the high-energy modes. 

We implement the above strategy in the Hubbard 
model on a 8 x 8 lattice, which has the set of one-particle 
states in momentum space represented in Fig. [j]. In 
order to study the effects of the enhanced scattering at 
the saddle-points, we start from a Fermi sea of occupied 
states below the shell e = 0, whose particle content cor- 
responds to a doping rate 5 m 0.22. We suppose that at 
a previous renormalization group stage the high-energy 
electron modes have been integrated out down to the re- 
gion between energies —A and A about the Fermi level 
shown in Fig. ^J. The modes integrated out have approx- 
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imately isotropic patterns about the top and the bottom 
of the band, and we may safely assume that their effect is 
encoded in a soft renormalization of the effective coupling 
for the low-energy modes, that we call U c s- Furthermore, 
we will deal with a situation in which U c s is in the weak- 
coupling regime, so that the low-energy physics can be 
described by building the many-body Hilbert space out 
of the one-particle states within —A and A. One has to 
bear in mind, however, that the interaction of the states 
at the Fermi level shell e — is always strong and essen- 
tially nonperturbative. In fact, we know from theoretical 
arguments that the main physical features have to come 
from the enhanced scattering with momentum transfer 
(-7T, 7r). This effect is properly captured in our construc- 
tion, given the high degeneracy of states at the e = 
shell. 

We proceed then to the numerical diagonalization in 
the many-body space of states of the low-energy theory, 
for different contents of particles at the Fermi level shell 
e = 0. We keep in all cases a Fermi sea of 24 particles 
filling the levels below e = 0. The instances with an even 
number of particles at the Fermi level are particularly 
interesting, since these are situations that allow to test 
the tendency to pairing in the model. Actually, there is 
a correspondence between the states with broken sym- 
metry in a macroscopic system and the sectors labelled 
by different quantum numbers in the finite cluster [ p"!] . 
The antiferromagnetic phase is characterized by a ground 
state in the sector A of spin singlets with total momen- 
tum P = (tt, 7t) and all the quantum numbers equal to 
I for the rest of discrete symmetries of the lattice. The 
d-wave superconducting phase may be characterized by 
a ground state in the sector D of states with total mo- 
mentum P = (0, 0) and odd parity under the exchange 
of the two axes in momentum space (the rest of quantum 
numbers remaining equal to I). A ground state belong- 
ing to the sector P of spin singlets with total momentum 
P = (0, 0) and even parity for the rest of the symmetry 
operations implies a phase without symmetry breaking 
at the macroscopic level. 

An important observation is that, for any even number 
of particles, the lowest-energy states found in the A, D 
and P sectors by considering the scattering of the one- 
particle states in the e = shell (and no particle-hole 
excitations from the Fermi sea) are always degenerate. 
Let us call = (7r,0), pb = (0, tt), pc = (37r/4,7r/4), 
Pd = (tt/4, 3tt/4). For two particles at the Fermi level, 
for instance, the following states have the same energy 

\A) = ±= ( C +(pa)c+( Pb ) - c+(p A )c+(p B )) \FS) 

~ iTf E (4(Pc)c+( Pc ) - e+(pc)4(p D )) \FS) 
\D) = ± ( c +(p A )c+( PA ) - c +(p B )c+(p B )) \FS) 



D ') = 7| E (4(pc)4(- Pc ) - c+( PD )c+(-p D )) \FS) 



n 



IP) = 



-j, (c+(pa)c+(pa) + c+(p b )c+(pb)) \FS) 
^ E (ct(PcH(-Pc) + c+(p D )c+(-p D )) \FS) 



R 



where 1Z stands for all possible reflections about the axes 
and \FS) represents the Fermi sea below the e = shell. 
It can be checked that this degeneracy among the lowest- 
energy states in the three different sectors holds for any 
even number of particles in the Fermi level e = 0, when 
particle-hole excitations from the Fermi sea are frozen 
out. 

Thus, in order to find to which sector belongs the 
ground state of the system, it suffices to make perturba- 
tions around the zeroth-order degenerate states. Wigncr 
theorem establishes that no generic level crossings are 
found when just one parameter is varied in a quantum 
system |p!q| . In our problem, we may switch adiabatically 
by means of an energy cutoff the number of particle-hole 
processes that are allowed in the space of states, in order 
to resolve the above degeneracy. The theorem guarantees 
that the ground state obtained in this evolution has to 
be adiabatically connected to the ground state of the full 
many-body space. In practice, it is enough to perturb 
the many-body states with up to two particle-hole exci- 
tations from the Fermi sea comprising 24 particles below 
e = 0. At this level, we may consider our computa- 
tional method as accurate as the coupled-cluster-double 
approximation, in which the relative errors estimated for 
the ground state energy of the Hubbard model are as 
small as ~ I0 -4 , for a 8 x 8 lattice and U = t |16|. As 
long as we are going to obtain our main results from the 
energy balance between different states in the model, we 
may expect them to be affected by smaller relative errors 
at that value of the interaction. 

We have computed the lowest-energy states in the sec- 
tors with different quantum numbers, within the above 
approximation to the space of states, for a number N of 
particles in the e — shell ranging from zero to 4. These 
values are about a doping rate S ~ 0.2. Otherwise, the 
case N = 4 already sets the computational limit in the 
number of particles that can be afforded with our Fermi 
sea. The ground state with all the closed shells at N = 
is a spin singlet with momentum P = (0, 0) and even 
parity for the rest of discrete symmetries. Its energy for 
U c tf = t turns out to be i=a — 25.I7I3i. The competition 
between the lowest energies in the different sectors for 
N = 2 and N = 4 is illustrated in Table |, for U eS = t 
and 2t. We have found that the ground state is always in 
the D sector of odd parity under the exchange of p x and 
p y . This conclusion for N — 2 can be considered a confir- 
mation of the significance of similar results obtained by 
the full numerical diagonalization of the Hubbard model 
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in a 4 x 4 cluster ]17| . In the present approach, the Fermi 
level at e — admits a larger number of particles, and 
it allows us to check the trend towards d-wave symmetry 
by filling the shell. 







P = (0,0) 
s-wave 


P = (0,0) 
d-wave 


P = (7T,7r) 


TV = 2 


(7 = £ 
(7 = 2t 


-24.7985 
-22.4508 


-24.7999 
-22.4540 


-24.7985 
-22.4507 


TV = 4 


U = i 
U = 2t 


-24.4270 
-21.7124 


-24.4273 
-21.7159 


-24.4253 
-21.7063 



TABLE I. Lowest energies of states with different quantum 
numbers, for particle content N = 2 and TV — 4 in the s = 
shell. 

To get insight about the nature of the ground states at 
TV = 2 and TV — 4, we have looked at the quasiparticle 
excitations supported by them. These correspond in our 
cluster to states with an odd number of particles (TV = 
1 and N = 3). States with N — 1 can be considered 
quasiparticle excitations over the fully symmetric state 
with TV — 0, and their dispersion may be obtained by 
looking at different momenta above the Fermi sea. We 
have plotted in Fig. || the lowest energies found for U c g = 
t among the many-body states with different momenta on 
the shell e = (including again up to two particle-hole 
excitations from the Fermi sea). 

The Fermi energy for N = 1 can be estimated by tak- 
ing the mean value between the ground state energies at 
TV = and N = 2, which turns out to be w -24.9856t for 
U c s — t. The lowest energy found at P = (ir, 0) is slightly 
above this value, by an amount of order 10~H. The 
quasiparticle excitations at P = f ) and P = (§, §) 
are found at higher energies. Results pointing in the 
same direction have been obtained in the diagonalization 
of small clusters |Q . We find a clear indication that the 
Fermi line is deformed by the interaction, in such a way 
that it is shifted inwards at the momenta (^p, f ) and 
(f;f) [@- This deformation is in agreement with the 
effect of pinning of the Fermi energy near a Van Hove sin- 
gularity that has been discussed on theoretical grounds 
pp|p|. According to this effect, the particles take ad- 
vantage of the strong screening of the interaction at the 
singularity to place themselves preferently near the sad- 
dle points of the dispersion relation, as observed in our 
computation. A similar phenomenon has been advocated 
recently in Ref. pjj . 

It is also of great interest the study of the lowest-energy 
states for TV = 3 particles on the e = shell. These 
may be considered as the quasiparticle excitations of the 
ground states for TV = 2 and TV = 4, and they are essen- 
tial to characterize the pairing of particles at even filling 
levels. The energies for the momenta on the e = shell, 
obtained in the same approximation as before, are plot- 
ted in Fig. | for U cS = t. 

The position of the Fermi energy at TV = 3 can be 




(tt,0) (3tt/4,tt/4) (tt/2,^/2) 

FIG. 2. Lowest energies of states with odd number of par- 
ticles N — 1 (squares) and TV = 3 (triangles). The values are 
given relative to the Fermi energy estimates in the text. The 
circles stand for the TV = 3 data after correction of the Fermi 
line shift. 



estimated from the ground state energies at TV = 2 and 
TV = 4, and it becomes w -24.6136i for U cS = t. The 
deviation from this value of the lowest-energy states for 
the odd-filled shell may signal the momenta at which a 
gap opens up for the quasiparticle excitations. Given 
the above mentioned shift of the Fermi line, though, one 
has to bear in mind that the Fermi level is not reached 
precisely at the points P = (^,f)andP = (|,f). This 
effect can be corrected by subtracting from the respective 
energies the shift from the Fermi level as measured for 
the TV = 1 quasiparticles. The final result can be seen 
m Fig. H The energy of the quasiparticles near P = 
(if ' f) ana ~ ^ = (§ ' § ) t urns out to be the same, within 
the precision of our computation, and slightly below the 
estimate from the mean value of the ground state energies 
for TV = 2 and TV = 4. On the contrary, the quasiparticle 
at P = (tt, 0) has an energy well above the Fermi energy 
estimate. It is clear that a gap opens up for quasiparticle 
excitations, but only at the position of the saddle points. 
The value A of the gap for U cS = t, A w 0.0017J, is 
sensibly greater than the energy difference between the 
lowest-energy states with s-wave and d-wetve symmetry 
for TV = 4. 

The conclusive evidence that the d-wave ground state 
at TV = 4 is built out of pairing of particles is given by the 
measure of pairing correlations with c?-wave symmetry. 
We have computed the correlations of the operator 
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5 (k + q)+ ff (k) 



3 +(k + q)c+(- 



k) 

(1) 



with g(p) = cos(p x ) — cos(p y ), in the ground state for 
TV = 4 at U c ff = t. The results are shown in Table 0. The 
momenta that can be measured are enough to discern a 
high peak at q = of the pairing correlations, together 
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with a sharp decrease away from the origin. To establish 
some comparison, we may quote that the peak of the s- 
wave pairing correlations in the lowest-energy state of the 
P sector is w 6.9. Taking into account that the d-wave 
correlations decay by one order of magnitude across two 
lattice units in momentum space, the translation of the 
results to real space is that the binding of particles is 
correlated on the whole dimension of the cluster. 
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q 


(OsCD(q)Oj CD (-q)> 


(0,0) 


25.969 


(f,0) 


4.907 


1 TV TV \ 

U' 4/ 

(f.o) 


8.526 


3.262 


/ TV TV \ 

\ 2 ' 4 ) 


1.671 


I TV TV \ 

\ 2 ' 2 ) 


1.059 



TABLE II. Values of the pairing correlator for different 
momenta near the origin. 



Although our results refer to a 8 x 8 cluster, they are 
robust enough to survive the extrapolation to larger lat- 
tices. The reason for that is the high degree of sym- 
metry of the model before switching on the interaction 
with the Fermi sea. We have seen that, by focusing on 
the scattering processes at the Fermi level shell e = 0, 
the lowest-energy state with momentum P = (tt, tt) and 
those with d-wave and s-wave symmetry turn out to be 
degenerate. This is a universal feature of the interaction 
of states on the nested Fermi line, irrespective of lattice 
size. Our approach to the many-body problem starts by 
solving first the strong interaction in the highly degener- 
ate Fermi level. While the interaction with the rest of ex- 
citations from the Fermi sea is relevant, as it resolves the 
mentioned degeneracy, its effect may be assessed within 
our weak-coupling renormalization group scheme. 

The existence of the d-wave ground state with remark- 
able pairing correlations provides a microscopic basis 
for a long-standing conjecture. There has been much 
hope that a mechanism based on antifcrromagnetic fluc- 
tuations may be the source of superconductivity in the 
copper-oxide materials. A different approach — though 
similar in essence — has put forward the idea that the 
Hubbard model should bear a form of Kohn-Luttinger 
superconductivity near half-filling. By stressing the rele- 
vance of the modes near the Fermi level shell at e = 0, our 
explicit construction of the d-wave ground state points at 
the processes with momentum transfer Q = (tt, tt) as the 
source of an attractive channel in the Hubbard model 
with bare repulsive interaction. Nesting of the saddle 
points at the Fermi level is at the origin of the relevance 
of such processes, which single out d-wave superconduc- 
tivity as the dominant instability of the system. 

I thank F. Guinea for useful comments on the 
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the Spanish Ministerio de Education y Cultura Grant 
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